Harmonic functions for a class of integro-differential operators. 

Mohammud Foondun 



Abstract 

We consider the operator C denned on C 2 (M. d ) functions by 

d „9 r, s d 



i,j=l 1 ■> i=l 1 

+ / [f(x + h)-f(x)-l m<1) h-Vf(x)]n(x,h)dh. 



lR d \{0} 

Under the assumption that the local part of the operator is uniformly elliptic and with suitable 
conditions on n(x,h), we establish a Harnack inequality for functions that are nonnegative in R rf 
and harmonic in a domain. We also show that the Harnack inequality can fail without suitable 
conditions on n(x, h). A regularity theorem for those nonnegative harmonic functions is also 
proved. 
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1 Introduction 



Researchers are increasingly using integro-differential operators (or equivalently, processes with 
jumps) to model problems from economics and the natural sciences. For instance, geometric 
Brownian motion is a standard model for a stock price. But this model is sometimes not satis- 
factory because it does not take into account sudden shifts of the stock price. To model this, one 
would like to use a process with some jumps, to represent the stock price. So understanding the 
properties of those operators is very important. 

The purpose of this paper is to consider functions that are harmonic with respect to the 
operator C, where 



is defined on C 2 (M. d ) functions. This is a typical example of a non-local operator, in the sense 
that the behavior of the harmonic function at a point depends on values of the function at points 
some distance away rather than just at nearby points. In probabilistic terms, the local part of C 
corresponds to the continuous part of the process while the non-local part controls the jumps of 
the process. The jump kernel n(x, h) represents the intensity of jumps from a point x to the point 
x + h and will be assumed to be nonnegative. 

We prove a Harnack inequality as well as a regularity theorem for harmonic functions with 
respect to the operator C without assuming any continuity of the coefficients ay, 6j and of the 
kernel n(x, h). We say that a function u is harmonic with respect to C in a domain D if Cu = 
in D; we give a precise definition in Section 2. Roughly speaking, the Harnack inequality states 
that the values of a non-negative harmonic function are comparable in a region. In other words, 
for all x and y lying away from the boundary of D, there exists a constant C not depending on u 
such that 



We show, with the aid of an example, that if n(x, h) does not satisfy some suitable conditions, then 
a Harnack inequality fails, while under mild conditions on n(x,h), a Harnack inequality holds. 

Since the fundamental work of Moser on Harnack inequalities for second order elliptic [20J 
and parabolicf?] partial differential equations with bounded and measurable coefficients, these 
inequalities have become increasingly important. Major contributions to this area have also been 
made by Krylov-Safonov[16j and Fabes-Stroock[10j. While there has been a lot of research on 
Harnack inequalities for functions that are harmonic with respect to differential operators, not 
much have been done for non-local operators. It is only recently that these results have been 
obtained for harmonic functions associated with purely non-local operators; see [6], [I] and [8]. 





u(x) < Cu(y). 
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The techniques we use to prove the Harnack inequality in this paper are similar to those in [6] and 
in [JJ, but have their roots in [16] where a non-divergence form elliptic operator was considered. 

As for the regularity theorem, we show that there exist a £ (0, 1) and a positive constant C 
not depending on u such that for all x and y lying away from the boundary of D, the following 
holds 

\u(x) - u(y)\ < C||u||oo|x - y\ a . 

Continuity estimates of the above type have a long history. Morreypjj] proved such an es- 
timate for second order elliptic partial differential operators in divergence form with bounded 
coefficients. His result, which was proved in two dimensions only, was independently extended to 
higher dimensions by DeGiorgi[llj and Nash [23]. Another proof was later given by Moser[20]. The 
corresponding result for operators in non-divergence form was established by Krylov-Safonov[16j. 
In [7] and [5], the authors considered purely non-local operators and proved a regularity theorem 
using probabilistic methods. It is also interesting to compare our result with the one obtained 
by Mikulevicius-Pragarauskas [22]. They considered a parabolic integro-differential operators and 
obtained a continuity estimate. However, their result, when specialized to the elliptic case, is a 
bit weaker than our regularity theorem. In that paper, the jump kernel n(x, h) satisfies a stronger 
condition than in our paper. Moreover, our techniques are different. 

Another paper which is related to our work here is that of Song-Vondracek [26] . Their result is 
a Harnack inequality for some discontinuous process. However, the jump kernel considered there 
is that of a a-stable process. Our result thus holds for a much wider class of processes. Related 
work also include a Harnack inequality for subordinate Brownian motion which has been obtained 
in [21]. 

The local part of our operator £ is of non-divergence form. In a forthcoming paper [9j, we 
consider an operator whose local part is of divergence form and whose jump kernel is symmetric. 
In that paper, the problem will be framed in terms of Dirichlet forms and a Harnack inequality 
together with a regularity theorem will be given. 

After stating the results in Section 2, we prove some preliminary estimates in Section 3. In 
Section 4, we prove a support theorem which is essential to our method. The proof of the Harnack 
inequality and regularity theorem are given in Section 5 and 6 respectively. In Section 7, we 
show that if the jump kernel n{x,h) does not satisfy some suitable conditions, then the Harnack 
inequality fails. 

2 Statement of results 

We begin this section with some notations and preliminaries. We use B(x,r) for the open ball 
of radius r with center x. We also use | • | for the Euclidean norm of points in R^, for the norm 
of vectors and for the norm of matrices. The letter c with subscripts will denote positive finite 
constants whose exact values are unimportant. The Lebesgue measure of a Borel set A will be 
denoted by \A\. 
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We consider the operator C denned by (1.1) and make the following assumptions: 

Assumption 2.1 We assume that the diffusion part of the operator is symmetric and uniformly 
elliptic and that the bis are uniformly bounded. In other words, there exist positive constants Ai 
and A.2 such that 

(a) the diffusion coefficients satisfy the following 

d 

AiM 2 < Yl V^A^Vo ^ Vlfl 2 . yeR d ,xe R d , 

(b) 

sup H&Hoo < A 2 . 

i 

We let A/"(Ai, A2) denote the set of operators of the form (1.1) satisfying Assumption 2.1. 
Besides nonnegativity, the following assumptions will also be imposed on n(x, h). 

Assumption 2.2 

(a) There exists a positive constant K such that 

[ (\h\ 2 A l)n(x, h)dh < K, Vx G R d . 

JR d 

(b) For any r G (0,1], any xq G R d , any x, y G B(xq,t/2) and z G B(xo,r) c , we have 
n(x,z — x) < k r n(y,z — y), where k r satisfies 1 < k r < kr~@ with k and (3 being posi- 
tive constants. 

n(x, h) can be thought of as the intensity of the number of jumps from x to x + h. n(x, z — x) 
thus represents the intensity of the number of jumps from x to z. So Assumption 2.2(b) says that 
the probability of jumping to a point z is comparable if x, y are relatively far from z but relatively 
close to each other. In Section 7, we show that such an assumption is needed for the Harnack 
inequality to hold. 

Since our method is probabilistic, we need to work with the Markov process associated with C. 
Let VL = D([0, 00)) denote the set of paths that are right continuous with left limits, endowed with 
the Skorokhod topology. Let X t (uo) = oj(t) for u G Q and Tt be the right continuous filtration 
generated by the process X. We say a strong Markov process (P x ,Xt) is associated with C if for 
each x, we have ¥ x (Xq = x) = 1 and for each x and for each u G C 2 that is bounded and with 
bounded first and second partial derivatives, u(X t ) — u(X ) — J Q * Cu(X s )ds is a local martingale 
under F x . This is equivalent to saying that P x solves the martingale problem for C started at x. 
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We assume that the martingale problem is well posed. In other words, we assume that the a^s 
and bis are continuous so that there exists a unique solution to the martingale problem. We make 
sure that none of our estimates are dependent on the modulus of continuity of the OjS and b, L s so 
that one can then use an approximation procedure to remove the continuity assumptions. 

For any Borel set A, let 

T A = inf{t : X t G A}, r A = M{t : X t £ A}, 

be the first hitting time and first exit time, respectively, of A. We say that the function u is 
harmonic in a domain D if u(Xt ATD ) is a P^-martingale for each x G D. If u satisfies some 
regularity conditions and Cu = in D, it is easy to see that u is harmonic in D. Since our 
operator contains a non-local part, our process will be have discontinuities. We write 

AV=limX„ AX t = X t -X t _. 

aft 

Our first result concerns the continuity of harmonic functions. Note that our hypotheses do 
not require Assumption 2.2(b) to hold. 



Theorem 2.3 Suppose Assumptions 2.1 and 2.2 (a) hold. Let zq G M. d and R G (0,1]. Suppose 
u is a function which is bounded in R rf and harmonic in B{zq,R) with respect to C. Then there 
exist a G (0, 1), C > depending only on the A/s and K such that 

(I ^ — y I \ a 
— ^ — J , x,y G B{zq,R/2). 

Our main result is the following Harnack inequality. 

Theorem 2.4 Suppose Assumptions 2.1 and 2.2 hold. Let zq G M. d and R G (0,1]. Suppose u is 
nonnegative and bounded on W 1 and harmonic in B(zq,R) with respect to C. Then there exists a 
positive constant C depending on the A/s, k, (3, R and K but not on zq, u, or ||u||oo such that 

u(x) < Cu(y), x,y£ B(z ,R/2). 

Remark 2.5 A chaining argument shows that both results above hold if R > 1 with C = C(R) 
depending on R. Theorem 2.3 does not hold for R > 1 with a constant which is independent of R. 

Remark 2.6 For the Harnack inequality, it is essential that u be nonnegative everywhere. Kass- 
mann 113}/ has shown that a Harnack inequality can fail for functions u that are harmonic with 
respect to symmetric stable processes of index a and where u fails to be nonnegative everywhere. 



5 



3 Some Estimates 



We start off this section with a proposition which allows us to assume A2 = when necessary. 
The proof is very similar to that of Theorem VI 1.2 in [2]. See also |25j . Define 

d 2 f(x) 



2 ^— ' dxidxj 
1,3=1 J 



+ / [f(x + h)-f(x)-l {mi) h-Vf(x)]n(x,h)dh. 



Proposition 3.1 Suppose C € jV(Ai,A2). If there exists a solution, say F, to the martingale 
problem for C started at x where C is defined as above, then there exists a solution F to the mar- 
tingale problem for C started at x. 



Proposition 3.2 There exist constants c\ and C2 not depending on xq such that if r < 1, then 
^ x ( T B(x ,r) — c i^) — t r ~ 2 f or x £ B{xo,r) and hence 



\ T B{x ,r) < c 2 r 2 ) < 



Proof. Let u be a nonnegative C 2 function that is equal to \x — xq\ 2 for \x — xq\ < |, which is 
equal to r 2 for \x — xq\ > r and such that its first and second derivatives are bounded by cr and c 
respectively. Since P x solves the martingale problem, we have 

E x u(X tATB{xo r) ) - u(x ) =E X Cu{X s )ds. (3.1) 

J 

Let us write the operator C as C = C c + Cd where 

d 2 u{x) , / du(x) 



, > 1 ^ — s - .O u(x) ^ — A , , 



2 ■4-', 13 dxidxj r-f * dxi 

i,3=l J 1=1 

and 

Cdu(x) = J [u(x + h) — u{x) — l(|ft|<i)/i • Vu(x)]n(x, h)dh. 

Since the first and second derivatives of u(x) are bounded, we have C c u(x) < C3 for x E B(xo,r) 
and hence 

C c u{X s )ds < c 3 t. (3.2) 
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Now let us look at Cdu{x) for x € B(xo,r) 



\C d u{x)\ 



< 



+ 



u{x + h) — u(x) — hh\<ih ■ Vu(x)]n(x, h)dh 
u(x + h) — u(x) — l(|h|<i)^- " Vu(x)]n(x, h)dh 



\h\<\ 



\h\>\ 
h + h- 



[u(x + h) — u(x)]n(x, h)dh 



By our assumptions and the fact that the second derivatives of u(x) are bounded, we get 



\h\<i 



[u(x + h) — u(x) — l(|/j,|<i)^- ' Vn(x)]n(x, h)dh 



< C3 / \h\ \\D u\\ oc n(x,h)dh < C4, 
!\h\<i 

u(x + h) — u(x)]n(x, h)dh 
n(x, h)dh < C5. 

\h\>l 



>\h\>l 
< 2||it|| 00 



Hence we have 



C d u(X s )ds <c 6 t. This together with fl3J]) and (J3T2]) yield 

E^(X tArB(xor) )<c 7 t, 



and so from r ¥ x (t b ^ xo ^ < t) < E, x u(X t /\ TB ^ XQ r) ), we get the first part of the proposition. The 
second part is obtained by choosing t = \r 2 . □ 

We have the following Levy system formula: 
Proposition 3.3 // A and B are disjoint Borel sets, then for each x, 



y2~L(x s -eA,x s <=B) ~ / / lA(X s )n(X s ,u - X s )duds 



(3.3) 



is a ¥ x -martingale. 

The proof is identical to that of the purely non-local operator and can be found in [6]. 
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Lemma 3.4 There exist c\ and C2 such that if r < ^, 
(a) K x TB( Xo ,r) — c iT 2 f or x S B(xo, r/2) and 
W E:Cr i?(xo,r) < c 2 r 2 for x G B(x Q , r). 

Proof. By Proposition 13. 2\ there exists C3 such that ^ >x { r rB( X0 ,r) — C3r 2 ) < ^(^(3.^/2) < C3r 2 ) < 
^. The first inequality follows by writing 

E a, r B (xo I r) > c 3 r 2 P :c (r i j (:C)r) > c 3 r 2 ), 

and using the above. Now let us look at the proof of the second inequality. For simplicity we 
assume that C G A/"(Ai,0). The general case follows by using Proposition 13.11 (and a change of 
measure argument). Since ¥ x solves the martingale problem, we have 

f tAr B(x ,r) 

E x u(X tATB(xo r) ) - u(x ) =E X £u(X s )ds. (3.4) 

J 

As before, let us write C = C c + Let us choose a bounded smooth function u(x) so that 
u(x) = \x — xq\ 2 for x G B(xq,2) and u{x) equals some constant greater than 4 outside the 
ball B(xq,A). Some calculus shows that Yli j=i ®ij u ( x ) = Yli=i duu(x) and is a constant for 
x G -B(xo,2). This and the uniform ellipticity of the local part of C implies that there exists a 
positive constant C4 such that C c u{X s ) > C4 whenever X s G B(xQ,r). 
To deal with the non-local part, we write 



Cdu(x) = / [u(x + h) — u{x) — h ■ Vu(x)]n(x, h)dh 
l\h\<i 

u(x + h) — u(x)]n(x, h)dh 

\h\>i' 
= h+h- 

Note that for \h\ < 1, we have x + h G B(xq,3/2) for x G i?(xo,r), so by convexity and the fact 
that n(x, h) > 0, we obtain 1\ > 0. As for the second term, we have 

I2 > / [\x + h — xq\ 2 — \x — xo\ 2 ]n(x, h)dh 

J\h\>l 

> 0. 

The facts that x G B(xq, r) and \h\ > 1 imply that x + h ^ B(xq, r) which means that the integrand 
is always non-negative. Combining the above, we have Cdu{X s ) > whenever X s G B(xq,v) and 
hence 



E x / Cu(X s )ds > C4E x {t A rB( x „,r))- 

JO 
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Now the left hand side of (|3.4p satisfies: 

E x u(X tATBM ) - u{x ) < r 2 . 

Combining the above and letting t — > oo, we get the second inequality. □ 

Corollary 3.5 For eachp > 1, there exists c± depending onp such thatforr < ^ and x € B(xo,r), 

E*(tL J < cir 2p . 

Proof. Note that 

K x (t bm ) > t¥ x (T BM > t). 

Letting t = 2r 2 and using Lemma 13.41 we obtain W(Tg/ X0)1 .) > 2r 2 ) < |. If Ot is the shift operator 
from the theory of Markov processes, then by the Markov property 

W x (T B ( X0 , r ) >{m + l)r 2 ) < F B (r B (x ,r) > ™r 2 , r B(*o,r) ° mri > r 2 ) 

= E x [F x mr* (r BM > r 2 ); t b{x ^ t) > mr 2 ] 

< ^(r B{x0!r) >mr 2 ). 

By induction F x {tb{x ,t) > mr 2 ) < (^) m . The required result then follows easily from this. □ 

The following result was first proved in the continuous case by Krylov in [15]. Since then, this 
inequality has been extended for diffusions with jumps. See for instance Theorem III15 of [17J. 
What follows is a consequence of Corollary 2 of [22J. 

Proposition 3.6 Consider C € JV(Ai, A2) and suppose that Assumption 2.2(a) holds. Let x € M. d 
and R 6 (0, 1]. If (F x , X t ) is a solution to the martingale problem associated with C G M(A\, A2), 
then for any bounded measurable function f , the following holds: 

E , j TB{ ^ R) f{Xs)ds < iV^II/H^ (3.5) 

JO 

where N depends on the A^s, and K. 

Remark 3.7 In the above, the constant N depends on the non-local part only through the constant 
K . It also does not depend on the radius R but the upper bound does depend on the radius as shown 
in (1331) . 
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The following gives a lower bound on the probability that our process hits a set A before 
leaving a ball. However, the set A should fill most of the ball. We will extend this result in the 
next section. 

Proposition 3.8 Let e > and r G (0, 1/2]. If x G B(x , §), A C B(xo,r) and \B(x ,r)-A\ < e, 
then ¥ x (Ta < T B r XQir \) > p(r, e), where there exists some eo such that p(r, e) > for < e < €q. 



Proof. From inequality f|3 . 5f) . we have 

rT B (x ,r) 



l A c{X s )ds 









E* / 




Jo 



l(B(x ,r)-A)(X s )ds 



< Nr\B(x ,r) - A\ 1/d 

< Nre 1/d . 



So we can write 



T B(x ,r) 



< 



{^ 2 B M ) 1/2 (V x (Ta < T BM )) 1/2 + Nre l l d . 



From Lemma [3. 41 and Corollary 13.51 we have E x 

T B( Xo ,r) > cir 2 and ^ x t b{xq r) < C2r 4 . So the above 

yields 

/ \ 2 

c\r 2 — Nre l l d \ 



» X (T A < r BM ) > 



The proposition is then proved with p(r,e) = ( cir 1 ) 2 . 



□ 



The following will be used only in the proof of the Harnack inequality. So far, this is the only 
place where we use Assumption 2.2(b). 

Proposition 3.9 Under Assumption 2.2, there exists a constant c\ which depends on K, such 
that if r < 1/2, z G B(xo, |) and H is a bounded non-negative function supported in B(xo,r) c , 
then 

^HiX } ) < Cl krE'H(X T ). (3.6) 



Proof. By linearity and a limit argument, it suffices to consider only H{x) = lc(%) for a set 
C contained in B(xo,r) c . From Assumption 2.2(b), we have n(w,v — w) < k r n(y,v — y) for all 
w,y £ B(xq, |) and v G B(xo,r) c . Hence we have, 



sup n(y,v — y)<k r inf n(y,v — y). 
yeB(x ,%) y£B{x ,%) 



(3.7) 
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By optional stopping and the Levy system formula, we have 



En ( Xt ^ S (, ,5) eC ) = E * h\x s -x s _\>%,x s ec)- 



E 2 / / n(X s ,v - X s )dvds. 

Jo Jc 



> E z (tAr B{x r } ) / inf n(j/,«-i/)A;. 

Letting i — > oo and using the dominated convergence theorem on the left and monotone conver- 
gence on the right, we obtain 

P^^GC)^^, I inf n(y,«-»)A;. 
Since EV^^r) > E z r b(«,|)i we have 

pZ (^ fl „ , 5 ) ^)> E*r B(-ii) / inf n(y,v - y)dv. (3.8) 

Similarly we have 

^ X °{X r } eC)< ^T B{x r } f SU p n(y, v - y)dv. (3.9) 

°' 5 JCyeB(x ^) 

Combining inequalities (|3.T[) . (|3.8|) and (|3.9|) and using Lemma [331 we get our result. □ 

Our process is a discontinuous one consisting of small jumps as well as big jumps. In many 
cases it is more convenient to discard the big jumps and add them later. This can be done by 
using a construction which is due to Meyer [18]. We will use this in the next section for the proof 
of the support theorem. 



Meyer's construction: 

Suppose that we have two jump kernels uq(x, h) and n(x, h) with no(x, h) < n(x, h) and such that 
for all x eR d , 

N(x) = / (n(x, h) — no(x, h))dh < c. 

JR d 

Let C and Cq be the operators corresponding to the kernels n(x, h) and n Q (x, h) respectively. If X° t 
is the process corresponding to the operator jCq, then we can construct a process Xt corresponding 
to the operator C as follows. Let S*i be an exponential random variable of parameter 1 independent 
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of Xt, let Ct = Jq N(X s )ds, and let U\ be the first time that Ct exceeds S\. At the time U\, 
we introduce a jump from Xu x - to y, where y is chosen at random according to the following 
distribution: 

This procedure is repeated using an independent exponential variable 52. Since N(x) is finite, 
this procedure adds only a finite number of big jumps on each finite time intervals. In [18], it is 
proved that the new process corresponds to the operator C. 



4 Support theorem 



The main result of this section is the support theorem. Before stating and proving this result, 
we present some ideas which will be crucial for its proof. More precisely, we will represent the 
solution of the martingale problem as a solution to a stochastic differential equation. We begin by 
representing the jumps of our discontinuous process as a function of a Poisson point process. 

Suppose that ¥ x is a solution to the martingale problem associated with C € Af(A\, A2) started 
at x. Let Y s be the point process associated with X s , that is, Y s = AX S if AX S ^ and 
otherwise. Then there exists a measurable function F(x,z) such that Y s = F(X S ^,Y S ) where Y s 
is a Poisson point process with intensity measure A. To make this statement more precise, let 
F(x, A) = {F(x, z) : z & A} and define 

Nt(A) = ^ l (AX s £F(X s -,A))- 
s<t 

Then under ¥ x , Nt(-) is a Poisson point process with intensity measure A. Moreover the measure 
A satisfies the following 

(a) 

(\z\ 2 A 1) X(dz) < 00, (4.1) 



l A (h)n(x,h)dh = I l A {F(x,z))X{dz). (4.2) 



/ 



Remark 4.1 The second condition above gives the relationship between the jump kernel n(x,h) 
and the Poisson process Nt(-). Moreover, the indicator function in \J^.^ can be replaced by a larger 
class of functions. For a more precise statement and proof of the above, see Theorem 12 in 
See also Chapter XIV of JT^. 
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We now relate (¥ X ,X), the solution of the martingale problem to that of a stochastic differential 
equation. Set fi([0, t] x A) = N t (A) and u([0, t] x A) = tX(A). Let W t be a Brownian motion with 
respect to the filtration J- t . Then X t solves the following stochastic differential equation 



dX t = <j{X t )dW t + b(X t )dt+ ( F(X t -, 

J\F(X t -,z)\<l 



z)(fj, — v)(dz, dt) 



+ [ F(X t _,z)v(dz,dt), X = x, (4.3) 

J\F(X t -,z)\>l 

where aa T has a™ as entries and a T denotes the transpose of a. The above has been taken from 
|17j . Chapter XIV of [12] contains more information about this relation. In fact, according to 
Theorem IIiq of [17], this representation holds under a more stringent condition on the big jumps 
of the process Yt (see Property M in [XTJ ) . Since the proof of the theorem below involves dealing 
with small jumps only and then adding the big jumps later, this does not affect our result (see the 
proof below). Here is our support theorem: 

Theorem 4.2 Suppose C G A/"(Ai,A2) and F x ° is a solution to the martingale problem for C 
started at xq. Let e > and suppose that (j) : [0,io] ~~ * ^> d differentiate with 0(0) = xo- There 
exist constants c\, C2 and C3 depending on A\, A2, to, K and sup i<io |</>'(^)l but not on e such that 
for all A > 0, 

r \ 2 t i 2 

¥ xo (sup \X t - <j){t) I < e) > a 1 - exp [ - Ae + — - (c 2 + c 3 e W )] . (4.4) 
t<t 12 J 

The support theorem says that the graph of X s stays inside an e-tube about (j). In other words, 
if = {(s,y) : \y — 0(s)| < e, s < t}, then {(s,X s ) : s < t} is contained in with a positive 
probability. 

Proof of theorem 14.21 We will assume C G 7V(Ai,0). The general case follows from Proposi- 
tion [3Tj We first consider the special case when the jump kernel of C is defined by no(x,h) = 
n(x, fo)l(|fti<i) and denote the corresponding process by X t . We will later use Meyer's construction 
to remove this restriction. We now use the stochastic differential equation representation of the 
solution to the martingale problem. In other words, we use the fact that Xt satisfies (I4.3p . Since 
our process do not have jumps of size greater than 1, the last term of (|4.3p can be taken to be 
identically zero. More precisely, replacing Ia(^) by | | Xf 1 /^i and n{x,h) by no(x,h) in (|4.2p . we 
obtain J^ F / X z - ) | >1 \F(x, z)\X{dz) = 0. Define a new measure Q by 

(h! ' exp [ - I ° ^'{s)o- l {X s _)dW s - \ f° ^'(s)^- 1 ^-)! 2 ^] . (4.5) 
1 Jo 1 Jo 

t 

J\F(X s -,z)\<l 



dF x o 

Let 



Zt = X t - f I _ F(X a -,z)(p-v)(dz,ds). 
Jo J\F(X,-,z)\<l 
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We see that 



J^<t>'{s)a- l {X s J)dW s ,z\ = /- j\'{s)a- 1 (X^)dW s , j^a{X s J)dW s 

cj)'{s)ds = -<j)(t) + 0(0). 



o 

So by Girsanov's theorem, under Q, each component of Zt is a semi-martingale. If 

W t = f a-^Xs^dXs- If a- x (X a J)F(X a -,z)(ji - u){dz,ds) 

JO JO J\F(X s -,z)\<l 

t 

-1/ 



a-\X s -)ct ) , (s)ds, 

then Wt is a continuous martingale and d(Wf, W^)=5ijdt under Q. Hence Wt is a d-dimensional 
Brownian motion under O. Note 



d(X t -<p{t)) = a(X t -)dW t + _ F(Xt-,z)(p-v)(dt,dz). (4.6) 

J\F(X t .,z)\<l 

If we prove the following: 

\X t -4>(t)\ <e) >C4, (4.7) 



t<t 

then the theorem will be proved, for if ^4 is the event {sup s<to \X S — 4>(s)\ < e}, then 

c 5 < Q(A) = [ (dQ/dF x °)dP x ° < (E X0 (dQ/dP a,0 ) 2 )3(P*°(i4))3. (4.8) 



The theorem then follows easily by noting the d<Q)/dF X0 has a finite second moment which is 
bounded by a constant depending on to, Ai and sup t<tQ |0'(t)[; see page 188 of [2]. Now let us 
look at the proof of (|4.7p . Let us write the left hand side of (|4.6p as dDt i.e, := — cp(t). 
Let A be a constant to be chosen later. Define 



\ 2 ft t*t r 

N t = XD t -— \a(X s -)\ 2 ds- I \ (e Az - 1 - Xz)n(X s ^, z)dzds. (4.9) 
2 Jo Jo J\z\<i 



i2 rt ft 

10 J\z\<l 

Set Y t x = e *. Then, by Ito's formula (for processes with jumps), we obtain 

y t A = 1 + X [ e Ns -dD t - f f e Ns ~(e Xz -1- Xz)n(X s ^,z)dzds 

Jo JO J\z\<\ 

+ J2[e N °-e N *--e N °AN s ). 

s<t 
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From Assumptions 2.1(a) and 2.2(a), there exist constants C5 and cq such that 

f \a(X s )\ 2 ds <c 6 t, 
Jo 

and 

\2 r 

(e Xz - 1 - Xz)n(X s ^,z)dz < — e |A| / \z\ 2 n(X s ^, z)dz 
\z\<\ 2 J\z\<l 

< otp*. 

By noting that AN S = XAD S , and using Theorem 10 of [T7] together with the above, we see that 
Y x is a martingale. The above bounds, together with (|4.9p also yield 

A 2 A 2 

Q(sup|A| > e) < Q(supe^ > exp[Ae - — c 6 t - — e l A lc 7 t ]). 

t<t t<t * * 

Since Y t x = e *, we can apply Doob's inequality as follows: 

A 2 A 2 

Q(sup |A| > e) < Q(su P y t A > exp[Ae - — c e t - — e W c 7 t ]) 
t<t t<t 1 1 

A 2 A 2 

< E Q Y t x exp[-Ae + —c 6 t + — e^c 7 t ] 

A 2 A 2 
= E Q Y X exp[-Ae + —c 6 t + — e |A| c 7 t ] 

A 2 A 2 

< exp[-Ae + y C6 t + Ye |A| c 7 t ]- (4.10) 

From the above we conclude that 

A 2 f 

Q(sup|A| < e) > 1-exp [-Ae + — ^(c 6 + c 7 el A l)]. (4.11) 
t<t * 

We now use Meyer's construction to recover the process X t so that (F x ,X t ) is a solution to 
the martingale problem associated with the operator C whose jump kernel satisfies the weaker 
Assumption 2.2(a). The trajectories of Xt now have jumps greater than 1. Recall that U\ is the 
first time that Ct exceeds S\ where S\ is an exponential random variable with parameter 1. More 
precisely, we have 

Q(sup \X t - 4>{t)\ < e) = Q(sup \X t - <f>(t)\ < e; U x < t ) + Q(sup \X t - <j>(t)\ < e; U x > t ) 

t<to t<t t<t 

> Q(sup \X t - <f>(t)\ < e)Q(Ui > t ). (4.12) 

t<t 
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Using the fact that 

Q(Ui < t ) < Q(5i < (sup AOi ) = 1 " e-( sup7V ) i0 , 

inequality (|4.12p reduces to Q(sup t<i() \X t — <ft(t)\ < e) > csQ(sup t<to |X t — 0(t)| < e) for some 
positive constant cs- This inequality, together with (|4.1ip and (|4.8p complete the proof. □ 



Remark 4.3 By taking A = e 2 in inequality 



we obtain upon choosing e small enough, 



P xo (sup \X t - 4>{t)\ < e) > ci[l - e- e3/2 ] 2 . 

t<to 

VFe now use i/ie /aci i/iai 1 — e _x > (1 — e~ l )x whenever < x < 1 to obtain 

¥ X0 { S u V \X t -(t,{t)\<e)>C2e\ 
t<t 



(4.13) 



for some positive constant ci not depending on e. However, C2 does depend on (ft via sup t<t \<p'{t)\; 
see (|4.8p and the discussion following it. 



We now present a corollary of the above support theorem. 

Corollary 4.4 Suppose C € AA(Ai,A2) and ¥ x ° is a solution to the martingale problem for C 
started at xq. Let e > and suppose that 4> '■ [0, to] ~~ > ^ *s continuous with 4>(0) = xq. There 
exist constants c\, C2 and C3 depending on e, A\, A2, to and the modulus of continuity of <f> such 
that for all A > ; 



r X 2 t i 2 

P IO (sup \X t - </>(t)\ < e) > ci 1 - exp [ - Ae + — -(c 2 + c 3 e |A| )] " 
t<t 12 



(4.14) 



Proof. Let us choose a differentiable function <f>& with derivative bounded by say C4 and such 
that sup s<t \4>(s) — (fid(s)\ < £• Moreover, we can choose 4>d(s) such that ||<^||oo depends only on 
t, e and the modulus of continuity of <p; see Page 60 of [TJ. Hence proving the following 



r X 2 t i 2 

P xo (sup \X t - <j> d (t)\ < e) > c 6 1 - exp [ - Ae + — ^(c 7 + c 8 el A l)] 

t<t L ^ 



will imply (|4.14p but with 2e instead of e. But the above inequality follows from Theorem 2? 
Hence the corollary is proved. □ 



Remark 4.5 The above corollary only requires the function (j) to be continuous but the downside 
of this generalization is that we can longer keep track of the dependence of the constants on e. 
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Let Q(x, r) denote the cube of side length r centered at x. If Ri denotes a cube with side length 
r, then Ri also denotes a cube with the same center but with side length r/3. The next result 
is not a probabilistic result. It enables us to decompose Q(0, 1) into smaller subcubes such that 
a subset A of Q(0, 1) fills a percentage of each of the smaller subcubes. Since this is Proposition 
V.7.2 of [2], we do not include a proof here. 

Proposition 4.6 Let q G (0, 1). If A C Q(0, 1) and \A\ < q, then there exists D such that (i) D is 
the union of cubes Ri such that the interiors of the Ri are pairwise disjoint, (ii) \A\ < q\Df)Q(0, 1)|, 
and (Hi) for each i, \AC\Ri\ > q\Ri\. 

A corollary of the support theorem is the following: 

Corollary 4.7 Let r G (0,R) and R G (0,1]. Let y G Q{0,R) with dist(y,dQ(0,R)) > r, 
C G jV(Ai, A2), and P be the solution to the martingale problem started at y. IfQ(z,r) C Q(0,R), 
then P(Tq( Z)7 .\ < Tq^ja) > C( r ) where £(r) > depends only on r, K and the AjS. 

The above two results together with the Proposition 13.81 are the main ingredients in obtaining the 
estimate below. The proof is essentially the same as that of Theorem V7.4 in [2] so we omit it 
here. 

Proposition 4.8 There exists a non- decreasing function ip : (0, 1) — > (0, 1) such that if B C 
Q(0, R), \B\ >0, Re (0, 1] and x G Q(0, R/2), then 

F*(T B <r Q{m )>iP(\B\/R d ). 

We now give a different version of the above proposition. This will allow us to use balls instead of 
cubes. 

Corollary 4.9 There exists a non- decreasing function (ft, such that if B C B(0,R), \B\ > 0, 
R G (0, 1] and x G 5(0, R/2), then 

F x (T B <T mR) )>cp(\B\/R d ). 

Proof. For simplicity, we assume d = 2. Higher dimensional cases differ only in notation. Let k 
be a large positive integer and let Rij be squares of the form [(£ — l)R/k, iR/k]x [(j — l)R/k, jR/k], 
where i, j G {— k + 1, — 1, 0, 1, k}. 

Take e > small enough and k sufficiently large so that 

C = {Rij : \Rij DB\>0, and Rij C B{0, (1 - e)R)} 

is nonempty. Let M be the number of elements in C. Let R*j be the cube with the same center 
as Rij but side length half as long. Let D = Up , e c5*,-. Pick z G R*,, where Rij satisfies 
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can choose k larger if necessary so that we can find such a cube. Then using 
Proposition 14.81 and the fact that \Rij\ = -p-, we have 

W z (T Rij nB<T Ri .) > ^j(\B\k 2 /MR 2 ) 

> H\b\/R 2 ), 

where the last inequality is obtained by noting that M < k 2 . Since dist(x, dB(0, R)) > R/2 and 
D C B(0,R), we can use Corollary 14.71 to obtain 

^(T D < r mR) ) > Cl , 

where c\ is a constant. Using the Markov property and the above inequalities, we obtain 

P x (T B <r mR) ) > W\P TR h(T B <r Rij );T Rh < r B{0>R) ] 
> Ci tP(\B\/R 2 ). 

□ 



5 The Regularity Theorem 

Now we are ready to prove the regularity theorem. 

Proof of Theorem 2.3. Let us suppose u is bounded by M in R d and z\ € B(z , R/2). Set 

r n = 2 p n , s n = 6ia n , for n G N, 

where a < 1, p < 1/2, and 6\ > 2M are constants to be chosen later. We choose 62 small enough 
that B(z 1 ,2r 1 ) C B(z ,R/2). Write B n = B(z 1 ,r n ) and r n = r Bn . Set 

M n = sup u(x), m n = inf u[x). 

x£B n x£B n 

We will use induction to show that M n —m n < s n for all n. The Holder continuity at z\ follows from 
this. Let tiq be a positive number to be chosen later. Suppose Mi — mi < Si for all i = 1, 2, n, 
where n > uq; we want to show 

M n+ \ — m n+ i < s n+ \. 

Let e > and choose z, y £ B n+ \ such that u(y) < m n+ \ + e and u(z) > M n+ \ — e. We will 
show that u(z) — u(y) < s n +i and since e > is arbritary, this will imply M n+ \ — m n+ i < s n+ i as 
desired. 
Let 

A n = {x £ B n : u{z) < (M n + m n )/2}. 
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We may suppose that |A n |/|S n | > 1/2, for if not, we can look at M n — u instead. Let A be a 
compact subset of A n such that |j4|/|i? n | > 1/3. Corollary 14.91 gives the following 

W X (T A <T n ) > Cl , (5.1) 

where c\ is a constant and x G B n+ \. Let z, y G B n+ \. By optional stopping, 

u{z)-u{y) = E z [u(X TA )-u(y);T A <T n } 

+ E z [u(X Tn ) - u(y);r n < T A ,X Tn G S n _J 
n-2 



+ ^ E*[tt(X T J - u(y); r n < T A , X Tn G B n ^ x - B n _ 



8=1 

+ E z [n(X T J - u(y);r n < T A ,X Tn £ S x ] 

= / 1+ / 2 + / 3 + / 4 . (5.2) 

By the Levy system formula, and Lemma 13.41 (see the proof of Proposition 3.5 of [B]) , there exist 
C2 and C3 such that 



sup ry(x Tn i B n _i) < sup w Tn [ 

eSn+i y&B n+1 J\h 



n(y, h) dh 

yeB n+ i y£B n+ i J\h\>r n -i-r n+ i 



= sup E y r n [ / n(y,h)dh+ / n(y,h)dh] 

y£B n +i J\h\>l Jl>\h\> 

The first term on the right of (|5.2[) is bounded as follows 

h < ( M "^ mn - J P y (TA < r n ) < ^ n P y (^ < r n ). (5.4) 

As for the second term, we have 

h < (M„_x - m^OF^ < T A ) < s n _i(l - V»(T A < r n )). (5.5) 



To bound the third term, we choose p = ^ A J i2Sc 3 ana - n °te that 



n-2 n-2 

a 

i=l i=l 

,2 



a 



a n {l - a)' 
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and 

P /a 



E2i ^ r P /" i 
Sn-i-lP <a n _i- 2 — 
i=i l-p 2 /a 



Using (|5.3p and the above, the third term is bounded by 

n-2 

i=i 

n-2 n-2 

< C 2 r^ ^ S n -i-l + C3 ^ S n -i-ip 2% 
i=l i=l 



a n (l-a) \-p 2 ja v 

128^ 



By our choice of p, we obtain 1 - p 2 /a > 3/4, p 2n /a n < l/2 2n and p 2 /a < jlft- so that the above 



reduces to 

n-2 

i=i 

a^ejci 4c 3 p 2 /a ] 
S Sn ~ l[ l-a + 3 J 
a 2 ela ci 

< Sn-lh + ™- 

1 — a 32 



We also choose 6*2 smaller if necessary so that # 2 < - 1 / — — 5 — — and obtain 

4 y 2a z C4 

ft < Sj ^- (5-6) 
Using (|5.3p again, we see that the fourth term is bounded by 

2M¥y(X Tn $ Si) < 2M[c 2 r2 + cgp 2 ^" 1 )] 
< 0i[c2a 4n ^ + C3 a 4n - 4 ]. 

By choosing no bigger if necessary and recalling that a < 1, we obtain for n > uq, 

h < ^=p- (5.7) 

Inequalities (|5.ip - (|5.7p give the following: 

u{y)-u{z) < \as n - 1 ¥ y {TA<T n ) + s n - 1 {l-¥ y {T A <T n )) + s n ^ 1 [^- + ( ^}. 

2 lb 8 
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Using the fact that a is less than one, we obtain 



u(z) - u(y) < 



_ Py(T A < T n ) Ci Ci" 

it, l 2 16 8. 



< Sn ri 5ci ' 
- a [ 16 



We now choose a as follows: 



This yields 



1 



5ci 
1~6~ 



u(z) - u(y) < s n a = s n+1 . 
The continuity estimate now follows from [20] . 



(5.8) 
□ 



6 Proof of the Harnack Inequality 

Proof of Theorem 2.4. By looking at u + e and letting e J, 0. We may suppose that u is 
bounded below by a positive constant. Also, by looking at au, for a suitable a, we may suppose 
that mfB( ZOl ij/2) it £ [1/4, 1]. We want to bound u above in B(zq, R/2) by a constant not depending 
on u. Our proof is by contradiction. 

Since u is continuous, we can choose Z\ G B(zq, R/2) such that u(-Zi) = |. Let rj = r\Ri~ 2 
where ri < ^ is a chosen constant so that Yli=i r « < -R/8. Recall that from Proposition 13.9^ there 
exists ci such that if r < |, y G B(x,r/4) and is a bounded non-negative function supported 
in -B(x, r) c , then 

^ X H{X TB(xr/2) ) < Cl k r EyH(X TB{xr/2) ). (6.1) 

For inequality (|6.ip to hold, we need Assumption 2.2(b). Let 7/ be a constant to be chosen later. 
Also, let £ be a constant defined as follows 

6 = - A — . 

5 2 Cl 

Let C2, C3 and C4 be positive constants to be chosen later. Once these constants have been chosen, 
we suppose that there exists x% G B(zq,R/2) with h{x\) = K\ for some K\ large enough so that 
the following is satisfied: 

— 1 — > 2, (6.2) 



k 



for all j. This is possible because of the fact that rj = r\Rj 2 . The constants k and are taken 
from Assumption 2.2(b). 
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We will show that there exists a sequence {(xj,Kj)} with Xj + \ E B(xj,rj) C B(xj,2rj) C 
B(z ,3R/4) with: 

= u(xj) and > i^e^'. (6.3) 

This would imply that Kj — > oo as j — > oo contradicting the fact that u is bounded. Suppose 
that we already have xi, x%, Xi such that (16.3H is satisfied. We will show that there exists 
Xi + \ € B(xi,ri) C 2rj) such that = u(x i+ i) and Kj+i > i^ie C2 ^ +1 ). Then by induction, 

(1631) will hold for all j. Define 

A = {ye B(xi,n/4)iu(3j) > ^^}- 

We are going to show that \A\ < — \B(xi, Tj/4)|. To prove this fact, we suppose the contrary. 

Choose a compact set A' C A with |.A'| > — \B(xi,Ti/4)\. Note that upon choosing n smaller if 
necessary, we can use (|4.13p to obtain 

^{T B{xhn/A) < r B(z0tR) ) > c 5 rl 

where C5 is independent of rj. To see this, consider (|4.13j) with e = rj/8 and let be a line 
segment joining z\ and x% (0(0) = zi and </>(io) = x«). Since zi,Xj € B(zq,R), \$(t)\ is bounded 
by a constant which is independent of i. Hence C5 is also independent of i; see (|4.8p and the 
discussion following it. 

Hence, using the strong Markov property, we can write 

F*i(T A , < T B{Z0)R) ) > E^[¥ Xt ^^(T a/ < T Bixijri) );T Bixuri/A) < r B(zoA) ] 
> <A((l/2) 2d+1 )c 5 rf. 

We now take C3 = 4> ((1/2) M+1 ) and C4 = C5. By optional stopping, the above inequality and 
the fact that u(X t /\T A , ) is right continuous, we obtain 

1 = u(zi) > E 2l [u(X rA , ATfl(zo ^);r A ' < T B{zQtR) ] 
> l^ P ^ (TA , <r ^ 0)Ji)) 
^ie C2 Vf +6 c 3 C4 
2. 
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This is a contradiction. Therefore \A\ < — \B(xi, r«/4)|. So we can find a compact set E such that 



E C B(xi,ri/4) — A and > ^\B(xi, r%/4)\. Let us write r Ti for TB( Xi ,r l /2)- From Corollary 14.91 
we have W Xi (T E < T ri ) > C6 where C6 is some positive constant. 
Let M = supg^. ,,^ u(x). We then have 

Ki = u(xi) = W. Xi [u(X TEATri );T E <r n ] 

+ E Xi [u{X TEATri );T E > r n , X Tri G Bfc, r,)] 

+ E x * [u(X TEATri y,T E > r ri , X Tri i B( Xi , n)] 

= I 1+ I 2 + I 3 . (6.4) 

Writing pi = ¥ Xi (T E < T r J, we see that the first two terms are easily bounded as follows: 

7i <i#|iii and h<M(i- Pl ). 

k 

To bound the third term, we prove E Xi [u(X Tr . ); X Tr . (jz B(xi,ri)] < r]Ki. If not, then by using 
(|6.ip . we will have, for all y € B(x il r i l '4), 



u(y) > Wu(X Tri ) >W[u{X Tr .);X TTi £ B(xi,n)] 

> ^ Xi HX Tri );X Tri * B( Xu n)] > ^ > 

contradicting the fact that < ^\B(xi, | . Hence 

h < r,Ki. 

So (|6.4f) becomes 



or 



M > 1 - 7/ - ^rf/fc 
Ki I- Pi 



1+ (i-frf/fc)»-i 

I- Pi 



Choosing rj = ^ and using the definition of £ together with the fact that pi > c% and rf? /k < 1, we 
see that there exists a positive L, such that inequality (|6.5[) reduces to M > -fQ(l + L). Therefore 
there exists Xj+i £ B(xi,ri) with it(xj + i) > + L). Setting -fQ+i = we see that 

= K-e log(1+i) . 
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The condition (|6.3p is thus satisfied provided we choose C2 = log(l + L). Finally, note that the 
fact that Yji=\ r i < f implies that B(xi, 2r { ) C B(z , 3R/4). □ 



7 An example 

In this section, we show that if an assumption along the lines of Assumption 2.2(b) does not hold, 
then the Harnack inequality can fail. This example is very similar to the one in [4]. But since 
we need some modifications and for the sake of completeness, we give a proof of the following 
proposition: 

Proposition 7.1 There exists a function n(x,h) which satisfies Assumptions 2.2(a) but not (b) 
and for which the Harnack inequality fails for functions harmonic with respect to the corresponding 
operator. 

Proof. Let B = 5(0,1), let y = (1/8,0) and for m > 4. let x m = (-1/8, 2~ m ), z m = (16,2" m ), 
C m = B(x m , 2~ m ~ 4 ), and E m = B(z m ,2~ m - A ). Define 

oo 

n(x,h) = l c m {x)l Em {x + h). 

m=4 

Note that n(x,h) satisfies Assumption 2.2(a) and not 2.2(b). Now we show that ¥ yo (Tc m < Tb) 
is small when m is large. We see that from Lemma 13.41 K yo Tg < c\ < oo. As before we are going 
to write C = C c + Cd- Now fix m, let e = 2 - " 1-4 , let g(x) = \x — x m \~^ ', where (3 € (0, 1) let 4> be 
a non-negative C°° function with support in 5(0,1/2) whose integral is 1, let <j) e (x) = e~ d (j){x/e) 
and let f £ = g* <j) e . Hence f e G C°° and we see that f e > C2e _/3 on C m . Since the local part is 
uniformly elliptic, we have \C c f e {x)\ < C3. From the definition of n(x,h), we have \Cdf e (x)\ < C4. 
Hence \Cf e (x)\ < C5. Since F yo is a solution to the martingale problem for £, 

E y °f e (X TCmATB )-f e (y ) = & f TCmATB Cf e (X s )ds 

Jo 

< c 6 K yo T B < 07. 

Hence 

c 2 e^F yo (T Cm < t b ) < E yo f e (X Tc . mATB ) < c 7 + f e (y ) < c 8 . 

Thus F yo (Tc m < tb) will be small if m is large. Now suppose that the Harnack inequality does 
hold for non-negative functions that are harmonic in B, that is, suppose there exists eg such that 

u{x)<c 9 u{y) x,y G 5(0,1/2), 
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for any nonnegative bounded function u which is harmonic in B. Let 

u m (x)=E x [l Em (X TB )]. 

Then u m (x) is bounded, nonnegative, and harmonic in B. Note that the only way that X Tg can 
be in E m is if X TB _ is in C m . We then have, using the assumption that the Harnack inequality 
holds, 

u m (y ) = W°[l Em {X TB )-T Cm <T B ] 

= E y " [E XT c m [l Em ( Xtb )];T Cm <r B ] 
= Ey°[ Urn (X Tc J;T Cm <T B ] 
< c 9 u m (x m )py°(T Cm <T B ). 



Then, 



(yo) c 9 pyo(T Cm <T B ) 

which can be made arbitrary large if we take m large enough. This is a contradiction and therefore 
the Harnack inequality cannot hold. □ 
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